Vacuum amplification 
of the high - frequency 
electromagnetic radiation 

G.A. Vilkovisky 

Lebedev Physics Institute and Research Center in Physics, 
Leninsky Prospect 53, Moscow 117924, Russia 



In memory of David Kirzhnits 



Abstract 

When an electrically charged source is capable of both emitting the electromag- 
netic waves and creating charged particles fi'om the vacuum, its radiation gets so 
much amplified that only the backreaction of the vacuum makes it finite. The re- 
leased energy and charge are calculated in the high-frequency approximation. The 
technique of expectation values is advanced and employed. 
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1 Introduction and summary 



Reaction of the vacuum on rapidly moving sources, or strongly variable fields is impor- 
tant for the evolution of black holes and early universe but is interesting also in electro- 
dynamics. We know that, in electrodynamics, the vacuum attenuates an external charge. 
Suppose now that the external source is not a monopole but, say, a dipole, and let this 
dipole be capable of emitting the electromagnetic waves so that the information about it 
reaches infinity. Then what will be the vacuum effect on such a dipole ? 

The answer obtained below is that the effect is opposite; a radiation of the dipole gets 
amplified. This effect becomes noticeable as soon as the typical frequency of the dipole 
exceeds the threshold of pair creation. A flux of charged particles that appears in this 
case is accompanied by an increase of the electromagnetic radiation. Generally, there is a 
nonlocal tail of radiation caused by the vacuum stress but, at a high frequency, the effect 
boils down to a multiplication of the classical radiation rate by a renormalization constant. 
Since the dipole is a nonlocal object, its renormalizationEI is finite and observable. 

The vacuum amplification of the electromagnetic waves emitted by a source is anal- 
ogous to the effect of the vacuum gravitational waves [1]. The difference is only in the 
theoretical mechanisms and in the dimensions of the coupling constants!. The dimension 
of the coupling constant causes that the gravitational effect never boils down to a mere 
renormalization. 

In the case of electromagnetic waves, the mechanism by which this effect emerges in 
theory is as follows. If one calculates the energy of charged particles created from the 
vacuum by a given nonstationary electromagnetic field [2], one finds that the result can 
be obtained only in the case where the electromagnetic field contains no outgoing waves. 
In the general case this energy is infrared divergent with the divergent term proportional 
to the energy of the outgoing waves. The appearing divergence is a signal that the 

^Residual after an infinite renormalization of the monopole. 

^The main difference is, of course, in the fact that, for a creation of the gravitational charge, there is 
no threshold. 
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calculation is not complete because the energy of the vacuum of charged particles goes 
partially into the coherent electromagnetic radiation. Indeed, the missing contribution 
comes from the backreaction of the vacuum on the electromagnetic field. If one calculates 
the effective electromagnetic field that solves the expectation-value equations, one finds 
that the quantum correction to the energy of the outgoing waves is also divergent, and the 
two divergences cancel each other. As a result, the total released energy is finite but is no 
more the energy of created charged particles alone. Rather it is a sum of the energy that 
goes with charged particles and the energy that goes with the enhanced electromagnetic 
waves. The two contributions can then be separated by calculating the released charge. 

It is worth noting that the vacuum reactions on the low-frequency and high-frequency 
external fields are very different. The effects like the anomalous magnetic moment in QED 
refer to the low- frequency electromagnetic fields and are not related to the effect considered 
here. On the other hand, in the mechanism described above one easily recognizes the 
physics that stands behind the so-called infrared disaster. Here this physics actually 
works, and, of course, there is no disaster if one considers not the transitions between 
concocted states but the evolution of expectation values. 

The terms " high-frequency approximation" and " strongly variable field" are used here 
as synonyms. Let I be the typical spatial size of the source of an external field and v 
be its typical frequency. On the other hand, let m be the mass of the hghtest particles 
interacting with this field. In the problem of the vacuum particle creation, the external 
field is considered as strongly variable if the energy hu dominates both the rest energy 
of the vacuum particle and its static energy in this field. The first of these conditions is 
discussed in Sec. 4 below, and its more accurate form is 



The second is exemplified in [2]. Under condition (1.1) the vacuum particles may be 
regarded as massless in the calculation of their fluxes. However, the mass m cannot be 
neglected in the calculation of the static polarization and charge renormalization. 

It makes sense to begin with quoting the result for the energy of particles created 
from the vacuum by strongly variable fields of arbitrary configurations. The respective 
calculation was carried out in [2] for the standard loop, i.e. for the vacuum action of the 




(1.1) 
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^,ac = ^logdeti^ , (1.2) 



form 



H = g'^'V^V,! +i^P- -Rlj - m'l (1.3) 

where the operator H is defined as acting on an arbitrary set of quantum fields. The hat 
over a symbol means that this symbol is a matrix in the space of field components, 1 is the 
unit matrix, and the matrix trace will be denoted tr. The external fields in (1.3) are the 
metric Qf^u, the matrix potential P, and an arbitrary connection defining the commutator 
curvature: 

[V^,Vj=7t^, . (1.4) 

The sign convention for the Ricci scalar R in (1.3) is such that, when acting on a scalar 
field, the operator H with P = and m = is conformal invariant. 

In the present paper, only the effect of the commutator curvature TZf^u is considered, 
and nevertheless the action (1.2) is needed with all the three types of external fields 
present. The dependence of ^vac on the metric is needed because the vacuum energy- 
momentum tensor is obtained by varying the action (1.2) with respect to the metric 

(and next using the retarded resolvent for the nonlocal form factors il). The dependence 
of ^vac on the potential is needed because the results for various quantum field models, 
e.g. for the spinor QED, are obtained by combining the standard loops with P generally 
depending on TZ^u (see [3] and Sec. 8 below). 

For the classical action of the commutator curvature one may take the expression 

^ci = J dxg^'hi n.A"" (1.6) 

with some coupling constant > 0. In the case of the electromagnetic connection, the 
is to be chosen so that (1.6) be the Maxwell action 



5', 



cl 



/ dxg'/'F^^F^'' (1.7) 



^See [2] and references therein. 
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(the matrix trace in (1.6) is always negative [2]). Denote 



The energy of the classical electromagnetic radiation and, in the high-frequency ap- 
proximation, also the outgoing flux of the vacuum energy can be calculated at the future 
null infinity X"*" [4,2]. The limit X"*" is defined as the limit of infinite luminosity distance r 
along the null geodesic that, when traced to the future, comes at the instant u of retarded 
time to the point of the celestial 2-sphere S. One has 



- -Ml^ci(.,0) + O(l) , (1.9) 
x+ du 



-;^^^vac(«,0) + O(l) (1.10) 



where 

Vv 



= Vu + 2Vr , {Vu, Vr) 

1+ 



= -l . (1.11) 



The notation d£ /du is introduced to represent the energy loss. Taken with the sign 
minus, each d£/du is the density of the respective outgoing flux of energy so that the total 
emitted energy is obtained by integrating {—d£/du) over the 2-sphere S (normalized to 
have the area 47r) and the time u. Specifically, the total released vacuum energy equals 

] duj d'S{ct>) =Y.eA (in vac |a+tall m vac) (1.12) 

— OO ^ / A 

and equals the total energy of particles created from the in-vacuum by external fields 
(see, e.g., [5]). Here a^^, a^^^ are the creation and annihilation operators for the out- 
vacuum, and Sa is the energy in the out-mode A. Since T^^ is energy-dominant [4], the 
flux {—d£c\/ du) is manifestly positive. The flux [—d£^/du) is sign-indeflnite because of 
the quantum uncertainty but the integrated flux (1.12) is positive [2]. 

Only the external fields generated by sources are considered in [2] and the present 
paper. The sources of external fields in (1.3) are 

J = P , Ji'^V^n^'' , J^'' ^R^"" -^g^"'R (1.13) 

where R*^^ is the Ricci tensor of the external metric, and the potential P is identified 
with its own source. These classical sources will be referred to as hare sources. The bare 
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sources are assumed to have their supports in a spacetime tube with compact spatial 
sections and a timehke boundary. Their domain of nonstationarity is assumed compact 
in both space and time [2]. 

At a large distance from a source, all its manifestations at both classical and quantum 
levels are governed by a single quantity, its radiation moment [2] defined as an integral of 
the source over a spacclikc hyperplane. The hypcrplane itself is defined as follows. One 
considers all timelike geodesies that, when traced to the future, reach infinity with one 
and the same value of energy per unit rest mass {E > 1) and at one and the same point 
of the celestial sphere {(f) e S). These geodesies make a 3-parameter congruence which 
is hypcrsurface-orthogonal, and the hyperplanes are the hypersurfaces orthogonal to this 
congruence [2]. Let 

= const. (1.14) 

be the equation of these hypersurfaces. The parameter 7 that, along with 0, labels the 
function T^^[x) is a redefined E: 

V£^2 - 1 

7= ^ , 0<7<1 (1.15) 
and the function T^^{x) itself is normalized by the condition 

(VT,^(x))^ = -(l-f) . (1.16) 
The radiation moments of the sources in (1.13) are the following integrals [2]: 

b = ^ j dxg^lH{T^^{x)-u)J{x) , (1.17) 
= j^Jdxg'/'6{T,^{x)-u)g"J^{x) , (1.18) 
D'^f' = ^ j dxg''^5{T,^{x)-u)g\g^-pr\x) (1.19) 

where g\ is the propagator of the geodetic parallel transport [6] connecting the integration 
point with the future end point of the geodesic having the parameters 7, 0. The moments 
are tensors at this end point depending parametrically on time u. 

At the limit 7 = 1 the hyperplane (1.14) becomes null. The vector and tensor mo- 
ments taken at 7 = 1 govern the classical electromagnetic and gravitational radiation. 
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Specifically, for the energy of the electromagnetic waves one has [2] 



d_ 
du 



£c\{u, (t>) = 



1 



d 



d 



7=1 



(1.20) 



The expansion of the vector and tensor moments at 7 = gives rise to the usual multipole 
moments [2] . The radiation moments integrated over 7 govern the energy of the vacuum 
particle production. One has [2] 

d 



+ Q.N. (1.21) 



where the abbreviation Q.N. means Quantum Noise and denotes the sign-indefinite con- 
tribution that is present in the vacuum energy flux because of the quantum uncertainty 
but sums to zero for the whole history [2,5]: 



J duj (fS{(t)) (Q.N.) = . 



(1.22) 



In the equations below the term Q.N. will often be omitted but its presence will be tacitly 
assumed in all expressions for the vacuum energy. 

Expression (1.21) is the starting point of the present work. It is seen that in the case 
of the vector moment (and only in this case) the validity of this expression is hmited by 
the condition 

(1.23) 



7=1 



which in view of (1.20) is a condition that the vector source does not radiate classically: 



d 

— £ci{u, 0) = . 



(1.24) 



If it does, the integral in (1.21) has a pole at 7 = 1. 

Below I shall consider only the contribution of the vector source assuming that the 
other contributions in (1.21) are absent. It is useful to decompose the vector moment 
over the vector basis at infinity: 
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(1.25) 



where m is the complex null vector tangent to the 2-sphere S, and m* is its complex 
conjugate. The projection 

(V^M + (1 - 7)V^r)z)'^ = e = const. (1.26) 
is the full conserved charge of the bare source [2]. Hence 



(1.27) 



The transverse projections of the moment taken at 7 = 1 define the complex news function 
of the electromagnetic wave£ 



ou 



____ = ICM) (1.28) 



SO that, by (1.20) and (1.27), 



Finally, the longitudinal projection of the moment 

V„rZ)" = Z)||(m,0) (1.30) 

plays no role in classical theory but, as shown below, it is responsible for the vacuum 
creation of charge. 

As pointed out in paper [2], the calculation in this paper is insufficient for a removal 
of the limitation (1.24) and needs to be revised. The revised calculation with all the 
needed amendments is carried out in the present paper. Its result is that the quantity 
(1.21) diverges indeed. The pole at 7 = 1 goes but its place is taken up by an infrared 
divergence. Only the total energy flux 

du du du 

is finite. The point here is that the electromagnetic field to be inserted in (1.31) should 
solve the expectation- value equations. To an appropriate order in the coupling constant. 



*I am using the terminology of Bondi [7]. 
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dSya.c/du can be calculated with the bare source but d£c\/du should already be quantum 
corrected, and this correction does not boil down to a renormalization of the coupling 
constant. This correction is infrared divergent, and it cancels the divergence in dSyg^/du. 
The final result is 



tot 



du 



1 1 ' 
(47r)2 3 



r 



u 

with the effective news function 
d 



tr gap 



7=1 



1.32) 



du 



Ces{u, 0) 



1 - 



,c - log 2 + 

24:71' ^12 



25.1 d 



du 



2A7rdu 



C{u, <P) 
d A 



J dr\og{m{u-r))—C{r,<t>) . (1.33) 



Here is the renormalized coupling constant, c is the Euler constant, and, as distinct 
from (1.21), the result is not independent of the mass m even in the high-frequency 

approximation. The quantities and C in the expressions above pertain to the original 
bare source. 

The retarded integral along X+ in (1.33) represents a nonlocal tail of the electromag- 
netic radiation. Technically, when a point tends to X+, its past light cone becomes a sum 
of a null hyperplane and a single null generator that merges with X+ [2]. The nonlocal 
radiation tail is a contribution of this generator (sec Appendix C). 

Since the energy flux in (1.32) belongs partially to charged particles and partially to 
the electromagnetic waves, the two contributions should, of course, be separable: 

dStot 



^^charge ^ ^^waves 



;i.34) 



du du du 

A calculation of the flux of charge helps to make this separation. The density of the 
outgoing flux of charge can be calculated from the expectation-value equations for the 
electromagnetic field (Sec. 3 below). The result for this density reads 



-|^e(«,0) 



K 



d 



3(47r)2 du 



7=1 



:i.35) 



so that the total released charge is 

I duj d'S{<P) [-— ej = / d'S{<P) = -oo) - = +oo)j 

(1.36) 



7=1 



Hence one may infer that the portion of the total energy flux (1.32) that goes with the 
charged particles is the one associated with the longitudinal projection of the moment: 



du (47r)2 3 



Jd^7'tr{-D,y . (1.37) 



The remaining energy in (1.32) goes with the electromagnetic radiation. Since anyway 
expression (1.32) is valid only in the high-frequency approximation, condition (1.1) can be 
used for its further simplification. It will be recalled that the domain of nonstationarity 
of the bare source is assumed compact. Its temporal scale (in time u) is a purely classical 
quantity of order l/i/. Therefore, if in (1.33) one writes 

log{m{u - r)) = log — + log(z/(M - r)) , (1.38) 

the contribution of the second term will be of order 0(1) whereas the contribution of the 
first term will be large: 

|.a,(„,,).(l-^,og^ + 0((^)°))|^C(„,« . (1.39) 

As a result, for u in the support of the bare news function, the radiation flux becomes 
merely a renormalized classical one: 

^%^ = Z^ , (1.40) 
ou ou 

Note the sign of the quantum correction! The radiation gets amplified. 

The results above pertain to the standard loop. For other models the vacuum fluxes 
are multiples of the respective fluxes for the standard loop (Sec. 8). Thus, for the spinor 

QED, the flux of charge is iAuice the one in (1.35), and the quantum correction to the flux 
of energy is twice the one in (1.32). Only the numerical constant which in (1.33) is 25/12 
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needs to be calculated anew but this constant is anyway unimportant. The explicit results 
for the spinor QED are obtained by introducing the said factor of 2 and substituting 

nf,, = -iqF^,i , e^-iqei , = V , tri = 4 (1.42) 

where F^i^ is the Maxwell tensor, e is the electric charge of the source, q is the electron's 
charge, and m in (1.33) is the electron's mass. 

In conclusion it will be noted that the result obtained cannot be the end of the story 
since, obviously, it violates the energy conservation law. Indeed, the frequency u is pro- 
portional to the energy of the bare source, and, since the factor logp/m can be however 
large, at a sufficiently large u the source will radiate more energy than it has initially. In 
this respect the present case is similar to the case of charged spherical shell considered in 
Ref. [2]. A spherical shell expanding in the self field emits no electromagnetic waves but, 
at an ultrarelativistic energy, its vacuum radiation violates the energy conservation law. 
The measure of the violation is in both cases one and the same, n^logu, and the cause 
is also one and the same: the problem has not been made fully self-consistent. Although 
in the present case the backreaction of the vacuum on the electromagnetic field is taken 
into account (otherwise the emitted energy would not even be finite), its reaction on the 
motion of the source is not. This task remains beyond the scope of the present work but 
it may be conjectured that the missing backreaction effect is nonanalytic in the coupling 
constant. 

Eqs. (1.32)-(1.36) and their corollaries are the main results of the present work. Their 
derivation is given below. A reader not interested in the technical details may still want 
to read Sees. 2 to 4. Sec. 2 presents the general scheme of the calculation including the 
important intermediate results and displays the mechanism of the vacuum backreaction. 
Sec. 3 presents the solution of the expectation- value equations and the calculation of the 
emission of charge. In Sec. 4, creation of massive particles is considered, and a criterion 
of the high-frequency approximation is derived. 

The technical details are presented in Sees. 5 to 8. The calculation required in the 
present work is more complicated than in [2] because the nonlocal form factors act now 
on functions having noncompact spatial supports. For a test function X, compactness of 
the spatial support is equivalent to the following powers of decrease at null infinities I"^ 
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and spatial infinity i : 




(1.43) 



The behaviours of the form factors derived or quoted in [2] are vahd only under conditions 
(1.43). For the present calculation, critical is the behaviour of the test function at X"*". 
The test function that doesn't satisfy condition (1.43) at X+ will be called singular at 
X+. The operators log(— □) and l/D with test functions singular at X+ are considered 
in Appendix C. The behaviours of the third-order form factors at are obtained in 
Appendix B. Appendix A summarizes the structure of the one-loop form factors. 
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2 The mechanism of the vacuum backreaction 



For obtaining the vacuum energy-momentum tensor to lowest order in the commuta- 
tor curvature, one needs the terms in the effective action quadratic in the commutator 
curvature and linear in the gravitational curvature, i.e. quadratic terms of order TZ x TZ 
and cubic terms of order R x TZ x TZ. Their general form is [8,9] 

'S'vac = 'S'vac(2) + S'vac(3) + highcr-ordcr terms , (2.1) 

^vac(2) = ^^^j dxg^'hi1Z,,^{-U)iz^^ , (2.2) 

'^vac(3) = j dxg^'hrY.^i{-n^,-n^,-U^)R{R^Mi) (2-3) 

with some form factors 7(— □) and rj(— Di, — ^2, — Da). In the basis of nonlocal invariants 
of third order [9] , there are 6 invariants of the needed type: 

R{R,2nM = Rr^TZ2a''TZ3p^ 

i?l7t27t3(18) = R^apyiAr^uilf 

R,n2nsm = Rr^^anr^VpTZs^, 

R{k2iz^{2Q) = R^v o^n^^y^n^p,. 

R{R2M'^i) = Rr^^Nxn^'^n^au (2.4) 

(I preserve the numbers that these invariants have in the full list of Ref. [9]). 

For in (1.6) to be the rcnormalizcd coupling constant, the form factor 7(— □) should 
satisfy the normalization condition 7(0) = 0. The normalized 7(— □) calculated for the 
standard loop is 

1 a) 
1 /" f 1—1 

^(-D) = -- da / dxlogil-x—) . (2.5) 



When applied to a high-frequency field, this operator takes the form 
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(The high-frequency hmit is considered in Sec. 4.) The constant 8/3 in (2.6) is observable 
since it accounts for the difference between the static regime in which the total initially 
stored charge is calculated and the high-frequency regime in which the emission of charge 
is calculated. Neither this constant nor the term in logm^ can be discarded when the 
operator (2.6) acts on a function singular at X"*" (cf. [2]). 

The third-order form factors Fj admit the massless limit and, in the high-frequency 
approximation, can be taken massless from the outset. One can then use the results of 
Ref. [9] where the massless Fj are calculated for all cubic invariants including the ones in 
(2.4). 

For obtaining T^^ at X+ one doesn't need the exact form factors. It suffices to have 
the asymptotic behaviours of Fj(— Di, — ^2, — Ds) with one of the arguments small and 
the others fixed. The difference with Ref. [2] is that these behaviours are now needed 
including the terms 0(n°). The algorithms of extracting the needed terms are derived in 
Appendix B. 

The contribution of the second-order action (2.2) to will be divided into two: 

= TZil) + (2.7) 
with T^{2) the contribution of the variation of the form factor: 

/ dxg'^XZi^) 6g,, -J^l dxg'/hr TZ,. 6^{-n) TZ^'^ . (2.8) 
Denoting T^^(3) the contribution of the third-order action (2.3), one has 

TZ-TZ{l)+TZm+TZ{S) . (2.9) 

The vacuum energy flux in (1.10) will then also be a sum of the respective three contri- 
butions: 

dSys^ ^ dSy^{l) ^ ggvac(2) ^ ^gvac(3) IQ) 

du du du du 

The expectation-value equations are obtained by varying the total action + Syg^ 

with respect to the connection field. These are the following equations for the source of 

the full, quantum-corrected, commutator curvature: 

JL^ + ^l{-0)J^n-J^^e (2.11) 
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with the retarded boundary conditions for 7(— □) [2]. Solving them iteratively one obtains 



■^FuU — "^Bare 27r^^ '-')'^Bare 



and hence 



'^FuU — '^Bare 27r '^Bare ■ 



(2.12) 



(2.13) 



For displaying the mechanism of the vacuum backreaction, it suffices to write down 
the expressions for T^^ and T^^(l) : 



rpflV 

-'cl — 



(2.14) 



r^(i) = --Ltr(7^'^S(-□)7^^-^^'^^7^<,^7(-□)7^"'') • (2.15) 



Using Eq. (2.13) one finds 



Pull 



•^="^Ba 



(2.16) 



The total energy-momentum tensor of the commutator curvature, T^j + , is then 



rpjJiV 

-'tot 



(rr + r^c(i) + rr(2) + 7r(3) 

(7;r-7rc(l) + TZ{2) + ITac(3)) 



Pull 



•^='^Ba 



(2.17) 



and hence the total energy flux is 



Q^tot _ /9£ci_9^vac(l) ^ d£y^{2) _^ d£y^{3) \ 

V rill, f)ii. rill. du ' 



'^='^Ba 



(2.18) 



du ^ du du du 

Thus the effect of the vacuum backreaction is changing the sign ofT^^l). As will be 
seen in a moment, this effect is dramatic. 

Note that if the substitution (2.13) was made in the action, then, after varying with 
respect to g^,^, both T^(l) and TJf^{2) would change their signs. This procedure is 
incorrect because it amounts to varying the action in g^y at fixed 7?.Bare whereas the 
energy-momentum tensor is obtained by varying the action in g^^ at fixed T^-puU- This 
makes difference since the relation between TZ^nW and T^Bare itself depends on the metric 
through the operator □. The correct procedure is making the substitution (2.13) in the 
energy-momentum tensor. 
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The dictum that at X"*" is infrared divergent means that expansion (1.10) does not 
hold. Rather there is an expansion of the form 



X+ 



logr 1 ^ / 1 \ 

terms — - — h terms + — ] 



(2.19) 



If this was the behaviour of the total energy-momentum tensor, the expectation-value 
spacetime would fail to be asymptotically flat. This is not the case but, in the inter- 
mediate expressions, the factor log r will conventionally be included in dE^^j du thereby 
considering this energy flux as divergent. 

The contributions (2.10) to dEvac/du are calculated in Sees. 5,6, and 7 below. Their 
main ingredient is the 7 = 1 radiation moment 



7=1 



(2.20) 



The latter notation is used everywhere below. The results are as follows: 



9^vac(l) 

du 



trj -(log mr + 2c - log 2 + 2 ) ) 



d 



(2.21) 



du 



—00 

-|-/*iog(„-x)(Ai,jM)(|i,,„M)} 



9gvac(3) 

du 



+ Q.N 
1 1 



(2.22) 



^ltr{-(log. + log2-^)(|-^?)(|-^,.) 



+ 1 J drlogiu-r){^Dt{r)){^DrM) 



d 



d 



2 



r 



1-72 



+ Q.N. . 



(2.23) 



Owing to the conservation law (1.26), only the transverse projections of Z)j survive in 
these expressions. Therefore, under the limitation (1.24) the contributions d8^{l)/du 
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and (9£vac(2)/9M vanish, and the contribution dSy3^{3)/du gives back the result of Ref. 
[2]. 

When the hmitation (1.24) does not hold, the contribution d£^{l)/du is infrared 
divergent. The contribution d£^{2)/du is not but it has another pathology. The total- 
derivative term in (2.22) does not vanish in the integral over time. On the contrary, the 
behaviour of this term at late time is 

— OO 

— oo 

so that the integrated flux (2.22) diverges: 

— oo 

The contribution d£^{?>)/du contains the divergences of both types. In the sum of 
the three contributions the divergence of the integral in time cancels but the infrared 
divergence doubles: 

Only in the total sum (2.18) with the changed sign of d£Yac{l)/du both divergences cancel, 
and the finite result (1.32) emerges. 

The cancellations outlined above do not depend on the relative sign and coefficient 
between and 5'vac (the in (1.6) is in fact kept arbitrary) but they depend crucially on 
the balance between 5'vac(2) and 5'vac(3). As seen from the expressions (2.21)-(2.23), there 
is a precise relation between the respective contributions, and this relation maintains 
for other field models (Sec. 8) despite the fact that S'vac(2) emerges from the purely 
electromagnetic coupling whereas S'vac(3) represents the vertices with the gravitational 
coupling. Owing to this relation, the final result is rigidly tied to the overall coefficient 
of the action 5'vac(2) which is merely the /3-function. A knowledge of this coefficient is in 
the end sufficient for obtaining the vacuum radiation fluxes. 
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3 The mean electromagnetic field. Emission 

of charge 



Since the quantum correction to the electromagnetic field cancels the infrared diver- 
gence in the vacuum energy, it should itself be infrared divergent. This point is clarified 
below but, before considering the expectation- value equations, it is useful to make a gen- 
eral analysis of the asymptotic properties of the commutator curvature and its source in 
the case where there is an emission of both waves and charge. To make difference with the 
notation already used, the quantities in this analysis will be distinguished with boldface. 

The existence of a fiux of charge at a large distance from the source implies that J" 
falls off at X+ like 

^y^^o(i,) (3.1) 



J" 



with some coefficient j°'{u, 0). This is the most general behaviour admissible for an isolated 
system. Although the support of the source J" is no more confined to a spacctimc tube, its 
domain of nonstationarity must remain compact in time in order that all fiuxes die out in 
the past and future of T+. More generally, the source should be asymptotically stationary 
in the past and future. To account for this property in the past, it will be assumed that 
the domain of nonstationarity of J" is confined to the interior of some future light cone 
u = U-. Then 

r{u,cj>) =0 . (3.2) 

U<U- 

The density of the fiux of charge from a source is expressed through the coefficient in 
(3.1) as follows [2]: 

1 



|;e(„.^) = ^V„.(r^j") 



V„Kj"(«,,^) (3.3) 



with Vv in (1.11). The function 

d 



e{u) = e(-oo) + J duj (fS{(f)) ^e(fZ, 0) (3.4) 
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defined by (3.3) up to an additive constant e(— oo) can be written as an integral over the 
future liglit cone [2]: 

e{u) = ^Jdx g^''^5{u - u)V-^u J'^(x) (3.5) 
provided tliat tlie constant 

e = e(— oo) (3.6) 
is taken as a conserved integral over an arbitrary (complete) spacelike hypersurface: 

e = / dxt'^6{T{x))VnT{x) J^(x) , (Vr)2 < . (3.7) 
47r J 

The function e.{u) may be called the Bondi charge, and the constant e the ADM charge 
since their meaning is the same as of the Bondi and ADM massed^. The ADM charge is 
the total charge of the source J in the initial state, i.e. before the beginning of emission. 
The Bondi charge is the charge that remains in a compact domain by the instant u of 
retarded time in the process of emission. The ADM charge is conserved because at any 
instant u it equals a sum of the charge emitted by this instant and the charge remaining 
by this instant, Eq. (3.4). 

Consider the conservation equation 

V„J" = . (3.8) 

Inserting the expansion (3.1) in (3.8) one obtains 

|^(V„Mr) = (3.9) 

whence, in view of (3.2), 

V„nr = . (3.10) 

The latter equation makes it possible to express the flux of charge in (3.3) through the 
longitudinal projection of j": 

-|-e(w,0) = i-V,rr(w,0) . (3.11) 

Thus the longitudinal projection V^r j" of the residue in (3.1) is responsible for the 
emission of charge, the projection Vq,u vanishes but no conclusion can be made on the 



continue using the terminology of the theory of asymptotically flat spaces [4]. 

19 



transverse projections j". Their vanishing does not follow and, at this stage, their role 
remains unclear. 

For obtaining the behaviour of the commutator curvature one must first consider the 
question of convergence of the moment Df of the source J. The analysis of convergence 
is carried out in [2]. When apphed to the present case, it gives the following result. The 
projection of Df on a basis vector in (1.25) converges if and only if the like projection of 
the residue j° vanishes. It follows that the projection V^m Df converges, the projection 
Dj diverges, and the behaviours of the transverse projections maDf remain unde- 
termined. Hence using the results for the retarded operator !/□ in Ref. [2] and Appendix 
C below one obtains 



= --Vaut)^{u,(l>) , (3.12) 

1+ r 



logr 1 



- J duVavr{u,(l>) (3.13) 



while for the transverse projections one has two cases: 



i) = , m«^J' 

ii) rriaT ^0 , rria^J' 



1+ r 

logr 1 



^-m„Df(ii,0) , (3.14) 



u 

j duma'f{u,(t>) . (3.15) 



T+ r 2 

— oo 

When the support of the source is confined to a spacetime tube, the convergent projec- 
tion VaiiDf equals the total charge of the source (Eq. (1.26)). This projection remains 
the conserved ADM charge also in the general case, even when there is an emission of 
charge and despite the fact that the integration hypersurface in Df is null: 

|-(v«KDf)=0 , V«uDJ = e . (3.16) 

The proof uses the explicit form of the null hypcrplane [2] and the stationarity of the 

source in the past. 

The fact that the projection VaV is generally divergent presents no real problem 
since this projection drops out of both the square of the differentiated moment in (1.27) 
and the commutator curvature. Indeed, solving the Jacobi identities to lowest order [2], 
one obtains for the commutator curvature 

Rm- = V.^J/, - V^^J. (3.17) 
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whence 



R 



(3.18) 

x+ ou ^ u ■ ■ u / \r^/ 
The projection (3.13) drops out of this expression by symmetry. Moreover, owing to the 
conservation law (3.16) one finds 



R 



2du 



{VuU - V^u my) (m* ^ J' 



(3.19) 



□ /J \Tf^ 

and the only projection of that can behave like 1/r (counting only powers) is 



-Vv m" R, 



T+ 



r du " (7-2) 



(3.20) 



where the coefficient dC/du at 1/r will conventionally be called news function although 
in the case ii) above it is infrared divergent. One has either 



i) rriar = , 



d 



d 



or 



fi) rn^'f 7^ , 



d - 1 
— C = (log r)-m„r («,(/)) 



(3.21) 



(3.22) 



At this stage there appears an argument to make a conclusion on the transverse fluxes 
m^j". If one wants the news function to be finite, these fiuxes must vanish. However, 
the only reason for insisting that the news function be finite is making finite the energy 
of the electromagnetic field since by (2.14) and (3.20) 



du 



IGttk^ 



tr(r V^v ruy R^"") (r V^t; R"^) 



(3.23) 



This is the reason indeed but only if the energy-momentum tensor of the electromagnetic 
field is given by expression (2.14). For a classical field it is. Therefore, for a classical 
field one has the case i) m^j" = 0, i.e. the only nonvanishing flux of J" is the flux of 
charge in (3.11), and the only divergent projection of the moment is VqV Df. The flux 
of the electromagnetic energy is then completely determined by the flnite news function 
in (3.21). However, if a c-number electromagnetic fleld is an expectation value rather 
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than the classical field, its energy-momentum tensor is not (2.14). Rather it is a sum 
rp^iu _|_ ^ g^j^^ ^j^g same argument that the energy should be finite may now be in favour 
of the case ii) j" 7^ where the news function is divergent. 

One is now ready to consider the expectation-value equations. In the high-frequency 
approximation, Eq. (2.12) takes the form 

^£^n = ^^a.e-^Q-log(-^)) J^3.e ■ (3-24) 

Since the bare source has a compact spatial support, one can use the result from Ref. [2]: 

2 d 



\og{-n)x 



-Di(«,0|X) + O(4) (3.25) 



1+ du ' ^r^ 

which is valid under conditions (1.43) and in which Di(-u, 0|X) is the 7 = 1 moment of 
the test source X. Since the local terms in (3.24) are 0(l/r^), one obtains 

1 d 



ja 
>^Pull 



x+ 127r du 



^1 Bare + 0(4) . (3.26) 



This is Eq. (3.1) with 

Hence, using the conservation law (3.16) for the bare source, one obtains 

^aUT = ^1 Baxe) = (3.28) 

and thereby checks condition (3.10). Next, using Eq. (3.11) one calculates the density of 
the flux of charge 

-li^.M) = -^,l{V^rDU..) (3.29) 

and thereby obtains the result (1.35). Finally, one calculates the transverse fluxes 

and discovers that they are proportional to the news function of the bare source. It follows 
that if the bare source emits waves, then the news function of the full source diverges. By 
(3.22), 

|;C«=(l-(log.)^)|;CB„ + 0(l) . (3.31) 
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However, one knows already that this divergence comes to the rescue. One can now check 
this again. From (3.23) and (3.31) one obtains 



cl 



du 



WV--^'°S'- + -=0(l))tr(|-CB„.)(|;C-U . (3.32) 



^FuU 

On the other hand, by (2.26), 

^-W<'°^^>'<S;^-=)(*^-)+0(l) . (3.33) 

As a result, the total energy flux (1.31) is finite. 

The approximate form (3.24) of the expectation- value equations corresponds to a 
neglect of the mass of the vacuum particles and is valid only in the region u > where 
the source is assumed strongly variable. This form can be used for a calculation of the 
derivative of the Bondi charge in the high-frequency approximation but cannot be used 
for a calculation of the ADM charge since the latter calculation involves the region u < 
where the source is static. The e/ecirostatic polarization with massless vacuum particles 
is infinite. Indeed, for the ADM charge (3.7) to converge, the full source must fall off at 
spatial infinity like 

j|iO=0(;^) (3.34) 
whereas a calculation with the massless form factor log(— □) in (3.24) yields the behaviour 



log(-n) Jeare 

and the divergent resultl 



o=0(i) (3.35) 



e Full = 1^1 - (log r)—j e Bare + 0(1) , m = . (3.36) 

The correct result for the ADM charge is obtained with the normalized massive form 
factor (2.5): 

CFuU = esare , m 7^ . (3.37) 



^The only exception is the case where the bare source has no monopole moment, e B^re ^ 0- Then 
one can show that also e -pnW — 0- An observable electric charge cannot be carried by massless particles. 
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4 Creation of massive particles and the 



high-frequency approximation 



In spite of their apparent similarity, the divergent renormahzation (3.36) of the ADM 
charge and the divergent renormahzation (3.31) of the news function have different status. 
The former is a result of an incorrect use of the high-frequency approximation in the static 
region whereas the latter is a natural consequence of the intensive pair creation. To show 
this and to derive a criterion of the high-frequency approximation, the expectation-value 
equations are considered below with the massive form factor 7(— □). 

The kernel of the operator (2.5) is obtained with the aid of its spectral form 



in which the threshold //^ = 4m^ appears explicitly, and, for the convergence at the upper 
hmit, the two spectral integrals are to be considered as a single integral. Outside the 
support of the source Jeare , Eq. (2.12) takes the form 



with the retarded resolvent (/i^ — O)^^. Here the order of integrations is important [10]. 
The spacetime integration implied in {/i^ — n)"^ Jsare is to be done first, and the spectral- 
mass integration next. 

One is presently interested in the behaviour of the full source at a large distance from 
the support tube of Jearc • At r ^ Z and u fixed, the retarded resolvent acting on a 
nonstationary source behaves as follows [10]: 



where f{u) is a positive function of time and angles having the dimensions l/u. 




(4.1) 




(4.2) 




(4.3) 
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1 °° 4^32 



oo 



Using (4.2) and (4.3) one can estimate the fluxes associated with created particles. 
For r ^ I one finds 

When projected on V^r, the coefiicient of in the latter expression is the density of 
the flux of charge through a tube of radius r. It follows that, because of the presence of 
the threshold, the flux through the tube of radius r » Hs suppressed by the factor 

Hence one infers that, although pair creation starts as soon as hu reaches the value of 
order mc^, the particles are created in the support of the source with small momenta and 
don't get far away. They stay in a compact spatial domain until hu reaches the value 

hu - mc' {^dj . (4.6) 

At this value there appears an observable flux of charged particles outside the support 
tube of J-Baie ■ The factor {mc/Ti)l may be interpreted as the number of created particles 
for which there is room in the spatial support of the source. If the creation is more violent, 
the particles get out of the tube. Finally, the high-frequency approximation is valid when 
hv is much bigger than (4.6): 

» mc^ f . (4.7) 



h 

Under this condition the flux of created energy and charge stops depending on the mass 
of the particles. The mass terms in (4.4) can then be discarded which is equivalent to 
replacing the second term of the spectral formula (4.1) with 

oo 

/V . (4.8) 

^ 

Hence the approximation (2.6) for the form factor. 

Expression (4.4) holds for the transverse projections of Jpuu as well: 



«4 ' (4-9) 
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and the coefficient of in this expression is nonvanishing whenever there is a nonvan- 
ishing flux of created particles outside the support of Jsare ■ Only in the special case where 
the electromagnetic radiation of Jeare is absent altogether may the transverse projections 
vanish. Thus the behaviour (4.9) is a direct consequence of pair creation. Then, by (3.22), 
the news function of the mean field inevitably diverges hke logr. Hence relation (3.31). 
The normalization scale of the logr in this relation can be read from the kernel of the 
operator !/□ (Appendix C). This is log(r/Z). 

It will be emphasized once again that there is nothing wrong about the mean elec- 
tromagnetic field. When pairs are created, its energy is no more governed by its news 
function since there appears a real vacuum contribution. The news function diverges but 
the vacuum energy redistributes and keeps the electromagnetic radiation down. 
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5 Calculation of T^^(l) at X+ 

vac \ / 



There remains to be presented the calculation of the energy fluxes (2.21), (2.22), and 
(2.23). This is, of course, the main part of the work. 

The T^^{1) is given by expression (2.15) with 7(— □) in (2.6). The commutator cur- 
vature to be inserted in T^^^ is the one generated by the bare source. The bare source has 
a compact spatial support. Therefore [2], 



□ 



(5.1) 



and hence by (3.18) 



r V Qu nil ' vr^ 



d_ 
du 



(5.2) 



T+ r 

It follows that for obtaining T^^{1) at one needs to know the behaviour of log(— □) 
with a test function that behaves at X"*" like 



X 



1+ r 

The needed result is obtained in Appendix C: 



(5.3) 



-iog( — Ax 



A{u,(j)) 



x+ 



(log mr + 2c — log 2) 



+ 



dT\og{m{u-T))^A{TA) + o(^) . (5.4) 



Substituting (5.2) for (5.3) one obtains 



(log mr + 2c -log 2) (|-D?)(|-£)i„) 



+ 



+ . (5.5) 



Here and below, use is to be made of the following identity: 



d 

du 



d 

j rfr log(M - r) /(r) = j cir log (m - r)—/(r) 



(5.6) 
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where /(r) is supposed to provide the convergence at the lower hmit. The convergence 
of the integral in (5.5) and similar integrals is provided by the assumption of asymptotic 
stationarity of the bare source. Under the simplified assumption that the domain of 
nonstationarity of the source is compact, there will be time instants w_ and u+ such 
that [2] 



d_ 
du 







u < 



d_ 
du 







U> Uj^ 



(5.7) 



In this way the result (2.21) is obtained. 
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6 Calculation of T^^(2) at J+ 



The T4'^(2) is defined by Eq. (2.8) with in (2.6). Using the spectral formula 

log(-^) ^- fd^i' - -^—^) (6.1) 

one finds 

oo 

(6.2) 

and the operator SO can be obtained by calculating 

I dxg^/HYn^P{5U)n^^ = I dxg^/^5g^,iY(^^^-{Vu^2)nfn2ap) (6.3) 

where p^'^(Vi, V2) is some polynomial in the derivative Vi acting on TZi and the derivative 
V2 acting on 7^2- In terms of this operator polynomial one obtains 

^^(2) = P'^^lVi, V2)i^§^^^7^r'^7^2a, (6.4) 

where the nonlocal form factor results from the spectral-mass integration in (6.2), and, 
up to higher orders in the curvature, the operators in (6.4) are commutative. 
The explicit form of Eq. (6.3) is 

j dxg^'Hxn^^{5u)n^p = I dxg^'^sg,, ix\^{y^n^^){y^n^p) - \g^''n{n''^n^p) 

-2Wa [n^pV'TZ'''^ - n^'^v'Tz^p) (6.5) 

whence 

TZ(^) - -77^2 fa ^"^^°^(°^^ VrA"^ V,^n2ap + total derivatives . (6.6) 
iz(47rj — U2 

A detailed analysis shows that the total derivatives in (6.6) either vanish at X+ or vanish 
in the integrated energy fiux, i.e. their contribution to d£^s^/du is quantum noise (Eq. 
(1.22)). The technique used in this analysis is outhned below. 
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In the remaining term of (6.6) use will be made of Eq. (3.17) and the conservation 
equation (3.8) to express the commutator curvature through its source: 



The first two terms of this expression are total derivatives. One obtains 
n^(2) = 



2a 



]_ J_ log(ai/a2) (^j derivatives 

6(47r)2 ^2 Di ^ 1 la 2 2 



(6.7) 



(6.8) 



where the indices /lu are symmetrized. Of the new total derivatives, the contribution 
of the first term in (6.7) is 0(l/r^) at X"*", and the contribution of the second term is 
quantum noise. The proof is given below. 

The form factor in (6.8) can be expressed through the operator Hq introduced in [2] 
and Appendix A below: 

i_ MEi^x.x. w . - / (l«,x,.,) ( / w) . (0.) 

—oo —oo 

The behaviour of this function as x — > T+ is obtained in the same way as in [2] by making 
the replacement of the integration variable 



q — r(T — u) , r — r{x) oo 



(6.10) 



where r is the new integration variable, and u = u{x) is the retarded time of the point x 
at T+. With q replaced as in (6.10), one has [2] 

HgXix) ^-Di{t,(P\X) (6.11) 
x^i+ r 

where the quantity on the right-hand side is the D-^ moment of the test source X. As a 
result, for (6.9) one obtains 



J_ log(Di/D2) 
□2 Dl - ^2 



XIX2 



^/,,(^Z>,(r.0|X,))(/^Z>x(..0|X.) 









f dr[ 


1+ J 





(6.12) 



Using the latter result in (6.8) one finds 



1 1 
—VuV'u — — rtr 



x+ 



6(47r)2 



r / 92 



dr d 



T — udr 

X) 

+ Q.N. 



DiJf] 



+0(1) . (6.13) 
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The integration by parts first in the internal integral 



I j^-Q=D^a{f)=\og{u-T)i^D^^{T)- / df log(M-f)^Di„(f) (6.14) 



d_ 



and next in the external integral yields finally the expression (2.22). 

It is now convenient to present a proof that the first two terms in (6.7) can be discarded. 
Their contributions to T^^^{2) are respectively 

1 



and 



J- vac 



'^-^ vac 



12(47r)' 



tr □ 



1 



6(47r) 



tr V«V^ 



^ log(e^ ^ ^ 



I logC./n.) ^,^^.^.^^, 



(6.15) 



(6.16) 



Using the same technique as above, one obtains 



1 log(ni/n2)^^^„^,,^^ 



1 



= VW^ii / 

2r 



df d 



T — M (?T 

— oo — oo 

Note that this behaviour is not even l/r^. It is 1/r. Nevertheless, one has 

°0(;)L^=0(1) (6.18) 

(see Appendix B). Therefore, the contribution (6.15) is indeed 0(l/r^) at X+. 

To calculate the contribution (6.16) at X"*", one may use the following result from [2]. 
If a symmetric tensor V^^^ is analytic at X"*", and its projections tangential to X+ vanish: 



Vo.uV'^^ =0(-V) , VauVpuV^^ =0(\v) , 



then 



(6.19) 



(6.20) 



By (3.16), the tensor (6.17) possesses the properties of V^^ . As a result, for the contri- 
bution (6.16) one finds 



^ vac 



x+ r 



(6.21) 
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where the respective energy flux is a total derivative in time: 



88. 



^ — 'vac 



-1 ^ / . 

tr — / ar 



du 12(47r)^ du J \J r — udr 



df d 



df d 
f — udf 



DUf)) . (6.22) 



oo — oo 



Since 



df d 



-oo — oo 



= -(Z)i(+oo)-Di(-oo))'^0 , 



(6.23) 



one obtains 



oo 

j du (z 



du 



. 



(6.24) 



Thus (6.21) is quantum noise indeed. 



32 



7 Calculation of T^^(3) at J+ 



The T^c(3) is obtained by varying the third-order action (2.3) with respect to the 
metric. Only the Ricci curvature that enters the basis invariants (2.4) needs to be varied. 
The commutator curvatures in these invariants are to be expressed through their sources 
via Eqs. (3.17) and (3.8). The result may be represented in the form 

rrc(3) = ^tr Y.'^i{-u,-n,,-u,)n,nr{i) (7.1) 

where the form factors r;(— □, —□2, — ^3) are linear combinations of rj(— □, —□2, — ^3), 
and the structures iZ2R'^^^{l) make some nonlocal tensor basis second-order in (the source 
of) the commutator curvature. The first argument □ of the form factors Ti is the oper- 
ator argument that in the action <S'vac(3) acts on the Ricci curvature. In the variational 
derivative (7.1) it becomes an overall operator acting at the observation point. In the 
diagrammatic language, the argument □ corresponds to the external line of the current 
(7.1). This is explained in more detail in Appendix B. 

It is shown in Appendix B that only the small- □ expansion ofri(— □,— □2,— 1^3) in 
the argument □ is relevant to the behaviour of the current (7.1) at X+. This expansion 
has the form 

f ,(-□, -02, -Ds) = log(-n)^K°2, Da) + Bi{n^, Dg) + 0{u) , □ ^ (7.2) 

and the contribution at X"*" of the terms 0(n) is already 0(l/r^). By the results in 
Appendix B, both Ai{02, O^) and B/(n2, Da) can be expressed through the operator Hq 
in a way similar to Eq. (6.9). Specifically, all Ai{02, ^3) are linear combinations of the 
following operators Fmn{^2, ^3)- 

^^,(□2, □3)^2X3 = -2 I rfgg™+"[(^)"'''7f,X2][(^)"^'7i,X3] . (7.3) 

—00 

Moreover, the term with log(— □) in (7.2) always appears in the combination 

log(-n)F„„(n2, 03) + L^n{02, Ds) (7.4) 
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with 



u 

j dq 



log(-|) + 2c)g-+- 



(7.5) 



All Bi{n2, Da) in (7.2) are linear combinations of Lmn{^2, ^3) and Fmn{^2, ^3). 
Thus one arrives at the ansatz 



nrc(3) 



^tr E(M-°H(°2,n3)+BKn2,n3))^2:^3''^(0 + C'(^) (7.6) 



T+ (47r)^ 

in which both the basis elements and the operator coefficients are presently to be de- 
termined. One can check that the following 9 structures quadratic in the source of the 
commutator curvature make a basis: 





— J2 ■ Ji" 


(7.7) 




= ^V'^V^(J2" • J3a) 


(7.8) 


n2nr{'i) 




(7.9) 


n2nr{^) 


= ^ V^V^V^V^ ■ J3^) 


(7.10) 


n2nr{^) 


= VaVM • V3'^V3V3'^) 


(7.11) 


n2V''{6) 


= V''{Vi'j2a-J3^) 


(7.12) 


7^27^3^"(7) 


= V''{J2a-Vtj,^^) 


(7.13) 


7^27^3'^"(8) 




(7.14) 


n2nr{'^) 




(7.15) 



The last two structures will be omitted since they cannot contribute to the energy flux 
through X"*". 

The respective coefficients Ai{p2-i 1^3) and Bi{p^2, 1^3) are obtained using the algorithms 
of Appendix B and the table of the third-order form factors in [9] . Only the basis element 
(7.10) with I — 4 has a nonvanishing Ai{^2, ^z)'- 



Ai{n2, 03) = , l^A 

^4(^2, 03) = -1 ^_(F22(n2,n3)-2Fn(n2,n3)) 



(7.16) 
(7.17) 
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This agrees with Ref. [2] where only the term with log(— □) in (7.6) was considered. The 
results for Bi{02, Da) are as follows: 



^1(02, Da) = , (7.18) 
^.(□.□3) = -i^ + l(i- + i-)FnC.n3) , (7.19) 

^3(02, Da) = l(i--i-](Fn(n2,n3)-Foo(n2,n3)) , (7.20) 



2 '-'3' 

i34(n2,n3) =i3'4(n2,n3)+i3"4(n2,n3) , (7.21) 

^'4(02, 03) = -i ^(L22(n2,n3)-2Ln(n2,n3)) , (7.22) 
B%{n,,o,) = l-L^F^,{^,,^,) , (7.23) 

O LJ2LJ3 



65(02, 03) = i-^(Fi2(n2,n3)-i^2i(n2,n3) 
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-Fn(n2,n3) -i^oo(n2,n3)) , (7.24) 
Be{02,a,) = _l(^i__3i-)Fn(n2,n3)-^ ^i^oo(n2,n3) , (7.25) 

67(02, 03) = ^^(i^ii(n2,n3)-^i^oo(n2,n3)) • (7.26) 

Only B4 contains the contribution of Lmn{^2, ^3) because the latter can appear only in 
the combination (7.4). 

The kernels for the superpositions of the operators Fjnn{^2: ^3) with 1/^2 and 1/^3 
are given in Eqs. (A. 15) and (A. 16) of Appendix A. By the same calculation as in Eqs. 
(6.9)-(6.12) their behaviours at X+ are expressed through the moments Di of the test 
functions: 



— -^mn('^2, n3)X2X3 
'-'3 



/ ,r(r-ur-'{^D,(rMX.)){^ J ^1X3)) . (7.27) 



Fmn{^2, n3)-^2-^3 



□2^3 



—00 —00 
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= 01X2)^1 (K^lXa) . 



(7.29) 



With these behaviours, the calculation of the Bi terms in (7.6) essentially repeats the 
calculation in Sec. 6. The Bi terms that involve the form factors (7.27) arc analogous to 
(6.8), and the Bi terms that involve the form factors (7.28) are analogous to (6.16). This 
concerns all Bi except B\. 

For I ^ 4 the results are as follows. Since Bi — Q, there remain two basis structures, 
/ = 6 and / = 7, in which the indices of the energy-momentum tensor do not belong to 
derivatives. Their contributions at X+ vanish by virtue of the conservation law (3.16). 
This may be exemplified with just one term 

^2 



1 



— Fn(n2,n3)J2aV3"J3^ 



□ 



- -;^V"« / ^r(r-«)(^Di.(r))(^Z)i(r)) 

—00 

+ 0(i)=0(i) ■ (7-30) 



The contributions of the remaining structures are of the form 



Bi{u2,n^)n2nr{i) 

with 



1 



1+ 



—c 



1 

=2~ 12 



_ 1 

=3~ 6 



_ 1 

l=5~ 4 



+ 0U (7.31) 



(7.32) 



The contribution of B'U is also of the form (7.31) with 0]^^^— 0. Thus the effect of all 
structures induced by the third- order action except the basis structure with I — 4 boils down 
to a finite renormalization of the classical news function. 

The main contribution comes from the basis structure (7.10) with / — 4. One may 
write 

log(-n)A(n2, a^)+BU{a2, a^ijAns'"' {4^) = V'^V'^(log(-n)/(x) +iV(x)) (7.33) 
where 



i{x) 

N(x) 



-l^aVp -^(^22(02, Da) - 2Fn(n2, , (7.34) 

4^"^^ 5^(^22(02, Da) - 2Ln(n2, . (7.35) 
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The scalar tr i{x) is the central object in Ref. [2]. In the same way as above one obtains 

d .X / d 



I{x) 
N{x) 



x+ 



x+ 



— oo 



(7.36) 



g^2 



+ 6^ / " O^S^^ - r) - log 2 + 2c + 1) (r)) (^^ia(r)) 

— oo 

+ + (7.37) 

where (7.36) reproduces the result in [2]. However, the scalar tr I{x) is needed at one 
more limit which in [2] is called i"*". This is the limit r — > oo along the timelike geodesic 
that, when traced to the future, reaches infinity at the point (f) of the celestial sphere with 
the energy E = [1 —7^)^^^^ per unit rest mass. The result obtained in [2] for this limit is 

7(1 — 7^) 7 , r d / d 



I{x) 



12r 



(7.38) 



where D" is the full (7-dependent) radiation moment of the source J. 

It follows from the properties of I{x) above that, for the calculation of (7.33) at X+, 
one needs to know the behaviour of log(— □) with a scalar test function that behaves at 
1+ like 



X 



1+ 



and at i+ like 



X 



7(1 - 7^) 



Q(7,</>) , Q(l,</')7^0 



(7.39) 



(7.40) 



where A{u, (j)) and (5(7, <l>) are some coefficients. The needed result is obtained in Ap- 
pendix C: 



-log(-n)X 



(7.41) 



and 



j cPS{4>)B{u,4>) 



u j d'S{<P) [-6g(l,0) +4 J d7^^(Q(7,0) -g(i, 

^ 

+ 0(logi^) . (7.42) 
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The behaviour of the function (7.33) at X+ is obtained by substituting (7.36) for (7.39) 
and using (7.41) and (7.37). Summarizing the calculation above one has 



rs;(3) 



= lv''''V"«(-%^)+o(i) , 



(7.43) 



^gvac(3) 

du 



{Any 



tr 



-^(logr + log2 + |)(Az,^)(p,„ 



d 



d 



B(m,0) + Q.N. 



The contribution of the latter total-derivative term to the radiation energy 

5£vac(3)- 



(7.44) 



(7.45) 



IS 



(47r)= 



tr 



d 



Substituting (7.38) for (7.40) and using (7.42) one obtains 



1 / d 



du d'S{(l)){-i—D'l)(—Di^ 



d 



—oo 

„2 



2^du ^^^du 



(7.46) 



(7.47) 



In this way the result (2.23) emerges. 
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8 Other models 



The results for quantum-field models other than the standard loop can be obtained 
by combining the results for the standard loop [3]. However, the results for the standard 
loop should then be known in full, including the contributions of the potential P. The 
contribution of the potential to the vacuum energy flux is given in expression (1.21) but 
this expression implies that the potential is regular at T+: 



0(3) • (8-1) 



x+ ^r" 

If condition (8.1) does not hold, the contribution of the potential should be calculated 
anew. For the starting point one may take the expression for 5'vac(2) with the form factors 
in the high-frequency approximation [2] 

^ ' 2{47r)^J ^ V 2 m^J 12 m^J 

+ const. PP + const. n^.^TZ'"' ] (8.2) 



and the expression for T^^^{3) calculated in [2] up to terms 0(n°) in the argument □ 
of the external line: 

TZ{^)\^ = (i^V'^V'^ tr log(-n)7(x) + 0(n°) , (8.3) 



i{x) = ^(Fn(n2,n3)-^i^oo(n2,n3))AA 

-\v^V^-^(F^2{^^,U,)-2F,,{n,^U,))j^-ji . (8.4) 

Of the missing terms 0{U^), the important ones can easily be restored. These are the 
terms in Lmn- Since the log(— □) in (8.3) originates from the expansion of the third-order 
form factors (Appendix B), each F^n in (8-4) should be accompanied by the respective 
L^n to form the combination (7.4). With the terms in added, expression (8.3) 
becomes analogous to (7.33). The remaining terms 0(n°) and the unspecified constants 
in (8.2) contribute only to a numerical renormalization of the news function (Sec. 7). 
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For the first example, consider the spinor QED. The effective action generated by the 
fermion loop in this model is (—1) times the action for the standard loop with 

p = \ri''n,. , YY + YY = '29'''i (8.5) 

and 

n^, = -iqF^A , tri = 4 (8.6) 

where F^^^, is the Maxwell tensor, and q is the electron's charge. 

Since the potential (8.5) is singular at X"*", one has to resort to Eqs. (8.2)-(8.4). One 
obtains 

ti P2P3 = -hi nrns^u (8.7) 

which is valid with any insertion of the form /(V2, V3). Then, by (6.7), 



tr P2P3 = tr 



1 / I c „, X /I c c flX 1 / 1 c i c „, 1 



When this expression is inserted in (8.4)-(8.3), the contribution of the first term vanishes 
at X"*" because of the presence of the overall □ (see Appendix B), and the contribution of 
the second term is pure quantum noise because it has the same structure as the J J term 
in (8.4) but with no form factor F22 0- As a result, one is left with 

i{x) = Fn(n2, ^.Wha - \v^Vp ^22(02, □3)^2" '^3'' + Q-N. (8.9) 

where the first term is the contribution of the potential. 

For I{x) in (8.9) the technique of Ref. [2] yields straight away 

/(x)L=-^ / (8-10) 

—00 

which is (—2) times the expression (7.36), and 



6r J ^du ^ ^du 



du{—D''){—D^) (8.11) 



-00 



^In each sum of Fnn in (8-4), only the Fnn with the highest n is to be retained since the junior Fn 
contribute only to the quantum noise [2]. 
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which is (—2) times the expression (7.38). Besides, there is the overall (—1) appropriate 
for fermions. It follows that dSyac{'i)/du for QED is twice the result for the standard loop. 
On the other hand, inserting (8.7) in (8.2) and changing the overall sign, one finds that 
5'vac(2) for QED is also twice the result for the standard loop. Thus, up to a numerical 
addition to the renormalization of the news function, all the results for QED are obtained 
by doubling the respective results for the standard loop and making the substitution 
(8.6). Note that, since the balance between 5'vac(2) and 5'vac(3) is maintained, the fact of 
doubling can be read just from the /3-function. 

Both the standard loop and the spinor QED have the "zero-charge" [11] sign of the 
static vacuum polarization. It is interesting to see what will be the results in the case of 
the "asymptotically free" sign. For that, consider creation of the Yang-Mills quanta in the 
external Yang-Mills field. In this consideration, it is convenient to refer to the standard 
loop with the commutator curvature 

n,^-K,u = c%F^. (8.12) 

where C|f, are the group structure constants, and is the strength of the external 
Yang-Mills field. 

In the minimal [3] gauge, the effective action generated by the ghost loop is (—2) times 
the action for the standard loop with P = and ^^i, in (8.12). 

The quantities pertaining to the loop of the gauge field will be distinguished with 
a tilde and expressed through (8.12). The loop of the gauge field is the standard loop 
with [3] 

P = pf^l^-^Ka.g'^ , (8.13) 

= '^%%,u = '^t,.€ ■ (8.14) 
Hence, in terms of the TZ^,, in (8.12), 

trPsPs = -AtrUrns^u , (8.15) 
tvTZrnsf^u = AtrUrn^^u . (8.16) 

Relation (8.15) differs from (8.7) only in the coefficient. Therefore, the calculation of the 
contribution of the potential repeats literally the one above; only the result should be 
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multiplied by 8. The contribution of the potential P to dSvac{3)/du is then (—24) times 
the result for the standard loop. The contribution of TZ/^i, is, by (8.16), 4 times the result 
for the standard loop. Since the contribution of ghosts is (—2) times the result for the 
standard loop, the grand total is (—22) times the result for the standard loop. The total 
action 5'vac(2) for the Yang-Mills field is also (—22) times the result for the standard loop 
as follows immediately from inserting (8.15) and (8.16) in (8.2) and adding the ghost 
contribution. 

Thus, also for the Yang-Mills couphng, all vacuum fluxes are multiples of the respective 
fluxes for the standard loop (with the substitution (8.12)), and the multiphcity is (—22) in 
accord with the /^-function but the price for the asymptotic freedom is that the radiation 
energy is negative. However, having taken off the head they don't weep for the hair [12]. 
Because the Yang- Mills quanta are exactly massless, a source of the Yang-Mills field would 
cause initially an infinite static polarization. The Yang-Mills charge is unobservable at 
infinity. 
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Appendix A. The one-loop form factors 



The basic building element for all one-loop form factors [9] is the operatoii 



l2q 



ng = \I^K^^^2qnj , q<0 (A.l) 

depending on the parameter q, with Ki the order-1 Macdonald function. By the properties 
of the Macdonald functions one has also 

d 



-n, = Ko [^2qn j (A.2) 
and 

2q—n, = an, . (A.3) 

Despite its scaring appearances the operator (A.l) has a simple kernel. Its retarded 
kernel is 

ngX{x) = — J dxg^/^6(^cx{x,x) -q)x{x) (A.4) 

past of X 

where a{x,x) is the world function [6], and the integration is over the past sheet of the 
hyperboloid of equal geodetic distance from x. The derivation of (A.4) is based on the 
spectral representation for the operator (A.l): 

n, = ---^q j ^/i ^ (A.5) 

^ 

where Ji is the Bessel function. Inserting in (A.5) the kernel of the retarded resolvent [10] 



1 If 1/0/ /iJi [n^/-2a , , _ 

X{x) = — / dxg^'H5{a)-e{-a)—^^j=^-L]x , (A.6) 



past of X 

doing the spectral-mass integrations 

J d^i Ji{ii^[^2q) = —1 



00 

J diinJi{ii\l ~2q)Ji{ii^ -2a) = 5{a - q) , (A. 7) 



*A11 operator functions are originally defined in the Euclidean domain □ < 0. 
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and using that 



(A.8) 



past of X 



one obtains (A. 4). 

The kernel (A. 4) was used in [2] without pointing out its relation to (A.l). This 
relation and the technique in [9] make it possible to obtain the kernels of all one-loop 
form factors in the expectation-value equations. Thus, for the second-order and third- 
order form factors [9] 



F (D a) - ( 9 Y ( d y togU,a,/na2] 



ii=i2=i 



Ja>0 



X 



(A.9) 



(A. 10) 



one has [9] 

rfemn(ni, ^2, Ds) 



{k + m + n)\ 

d^ 



dq q 



k+m+n 



X (y^) —^K, (v^) —Ko (v^) (A.12) 



and, therefore, 

F^^{o^,n2)XiX2{x) 



-2 / dqq-^-[{fy^\x,{x)][(fy^\x2{x)] , (A.13) 



'dq 



dq'' 



rw(ni,n2, □3)^1^2^3(2:) 



(A; -|- m -|- n)! 



dq q 



k+m+n 



X 



{^y'\x,i4 l{^r\x2ix)] [{^y'\x,ix)] (A.14) 

with HgXix) in (A.4). 

Eq. (A. 3) makes it possible to obtain easily the superpositions of the kernels above 
with !/□. For example, 

i-F..c., □2)x,X2 ^ - / dqq--\^,n,x,) (^^ / fn,X2 



(A.15) 
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□ 



which is vahd including the cases n = 0, m = (cf. [2]). Similarly for the third-order 
form factors. The convergence of the integrals in q at the upper limit is controlled by the 
behaviours 



He 



q-*0 



O(logg) 



(A.17) 



q=o n ' dq 

following from (A.l), and the convergence at the lower limit should be provided by the 
properties of the test functions [2]. Eq. (A. 3) can also be obtained directly by acting with 
the operator □ on (A. 4) and neglecting the curvature in a, Da = A + 0[R]. 
One might have introduced the kernel even more elementary than (A. 4): 



1 r 

egX{x) = ^ y f^^O{q - a{x, x))X\ 



x) , q < 



past of X 



d_ 

dq 



— Qq — Hq 



Or 



q=0 



□ 2 



whence 



e, = |x,(y2^H) . 

The initial condition in (A. 19) implies that 

r^(x) = ^ / dxf/'e{-a)X 

past of X 

and 

^ J dxf'H(a{x,x)) j d^f'H(a{x,^))X 



1 

□ 2' 



(A. 18) 
(A.19) 

(A.20) 
(A.21) 



past of X 



past of X 



1 f 

= — j dxf'^e(-a{x,x))X . (A.22) 

past of X 

Eq. (A.21) can be obtained by acting with the operator □ on (A. 18) and using (A. 8). It 
is also a limiting case of the formula [2] 



1 



{w? - □) 



^X{x) 



1 

Stt 



J dxg^/^e{-a)Jo(rnV^)x (A.23) 



past of X 
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for the massive operator. The kernels (A. 18) and (A. 21) do not decrease at the future 
infinity and can be used for a direct determination of the moments. Thus, 



du 



-Di{u,<P\X) 



(A.24) 



where the quantity on the right-hand side is the Di moment of the source X. 
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Appendix B. The third-order form factors at X+ 



The form factors in the third-order action (2.3) are hnear combinations of the functions 
Tkmni^i, ^2, 1^3) introduced in Appendix A. The typical contribution of such a form factor 
to the energy-momentum tensor at point x has the form 

rw(n,n2,n3)X2X3(x) (B.i) 

where the X's arc the commutator curvatures or their derivatives, and it is assumed that 
first ^2 acts on X2 — X{x2), and ^3 on X3 — X{x3) with subsequently making the points 
X2 and X3 coincident with the observation point x, and next the first argument □ of the 
form factor acts on the thus obtained function of the observation point. This nonlocal 
structure corresponds to the diagram in Fig. 1. 

By the results in Appendix A, expression (B.I) can be represented as follows: 

r.^.(n,n2, □3)^2X3(0;) = ^^^^^ J dq U'^,K, (^jj^(g,a;) (B.2) 

where the operator □ acts to the right on the function of x, and this function is 



/ d \m+l 



d \ n+1. 



^(g,x)=?"^+"[(-J n,X2(x)\[[-} n,Xs(x)\ . (B.3) 

When the X's are expressed like in Eq. (6.7) through the sources J having compact 
spatial supports, there occur two essentially different cases. An example of the first case 
is 

X2Xs = J^jf , (B.4) 
and examples of the second case are 

X2X3 = ^ J2"J3'' or J2"^^J3^ or V„V;3(^J2"J30 • (^•^) 

The difference between the two cases is in the behaviours of integrals with the function 
T{q,x) as X ^ J"*". These behaviours are readily obtained by the technique in Ref. [2] 
(see also Sec. 6 above). In the first case one has 





/ 



-00 



dqJ^{q, x) 



= 0{r-^) , (B.6) 
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and in the second case 



J dqj^{q,x) 



0{r-^) 



(B.7) 



The second case is our main concern here since the function (B.7) is singular at X"*". 

Our present goal is obtaining the behaviour of the current (B.2) as x — > T+. The 
principal assertion is that this behaviour is determined by the first few terms of the 
small- □ expansion of the form factor F in the argument □ of the external line. For the 
proof it suffices to consider two generic terms of the small- □ expansion of the function 
Ko {^/2^) in (B.2): 

{qn)P and {qa)P logqn . (B.8) 

It will be recalled that the behaviours (B.6) and (B.7) are obtained by making the 
replacement (6.10) of the integration variable q. From the form of this replacement, it 
follows that in the case (B.7) one has also 





/ dqqP^iq, 



X 



dqqP \og{—q)J^{q,x) 



(B.9) 

:b.io) 



and then, by the result in Appendix C, 
log(-n) ( J dqqPJ^{q,x) 



r-f-Mogr ap(«,0) +0(rP-2) , p>l (B.ll) 



with one and the same coefficient ap{u, 0) in all the expressions (B.9) to (B.ll). 
Moreover, using the following form of the operator □ at X+ [2]: 



DX 



1+ 



one obtains 



UO{rP) 

and, in the exceptional case p — —1, 



(B.12) 



(B.13) 



1+ 



0{r-^) 



(B.14) 
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Owing to the latter fact, one has 



p>i . 



:b.i5) 



The relations above make it possible to obtain the contributions of the expansion 
terms (B.8) to (B.2): 



r 1 

aP / dqqPT{q,x) = 0(—) , p>l 



u u 
log(-n)(^y" dqqPJ^{q,x)^ + J dq qP \og{-q)J^{q, x) 



(B.16) 



O 



P>1. 

(B.17) 

It follows that, in the case (B.7), the function of □ in (B.2) can be truncated as follows: 



d^ 

^'^^0 ( x/2gn 



-{-lf{k-l)\ , k>Q 

1 qU 
— log c , /c = 



+ irrelevant terms (B.18) 

where the irrelevant terms are the terms whose contributions to (B.2) are 0(l/r^) at . 
By a similar analysis, in the case (B.6) this function can be truncated even more: 

, fc> 

= 



2gn 



1 



2M-°) 

+ irrelevant terms 



:b.i9) 



and one recovers the algorithm used in [2]. Thus the amendment needed in the case (B.7) 
as compared to (B.6) is retaining the terms 0(n°) of the form factors. 

Besides the contributions of the form (B.l), the vacuum energy- momentum tensor 
contains contributions in which the form factors VkmniP) ^2, Da) are superposed with 
!/□ in the argument of the external line [9]: 



- r 

□ 



kmn 



;□, ^2, □3)X2X3(X) 



:B.20) 



These contributions occur only at /c > 1 [9] and only in the case (B.6) i. By the same 

consideration as above, the operator function in (B.2) can then be truncated as follows: 

^This fact is a matter of a direct calculation [9] but it is also a necessary condition for the expectation- 
value spacetime to be asymptotically flat. 
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1 t <'" 1, 



2qU 



1 



+ 



|(log^+2c-l)g 



k > 1 



+ irrelevant terms 



:b.2i) 



The effect of these truncations is that the third-order form factors boil down to the 

second-order form factors. The latter are the functions Fmn and Lmn introduced in Sec. 7, 
and similar functions originating from expansion (B.21) and differing from Fmn and Lmn 
by an extra power of q: 





Gmn{02,as)X2Xs^-2 J dq Q 



^dq^ 



[B.22) 



u 

[n^, a,)X2Xs = -2 J dq (log(-|) + 2c) q 



m+n+l 



'dq^ 



3 ■ 
(B.23) 



Using Eq. (A. 3), the latter functions can be expressed through Fmn and Lmn' 

2 

Gmn(02,03) = —{Fm+2A^2,Os) + (m + l)Fm+l,n(^2,Os)) 

= p^(i^m,n+2(n2,n3) + (n+l)F^,n+l(n2,n3)) , (B.24) 



2 

Mmn{02,Ds) = —{Lm+2,n{^2:03) + {m+l)Lm+l,n{02:03)) 

= ^(L„,n+2(n2,n3) + (n + i)L^,n+i(n2,n3)) ■ (B.25) 



Equivalence of the two forms in (B.24) follows from the identities for Fmn in Ref. [2]. 
Similar identities can be derived for Lmn, G^n and M^„. All of them are based on Eq. 
(A. 3) and the integration by parts in the q integrals. 

The consideration above can be summarized as follows. For the case (B.4) one has 



^kmn{0,02,03)J2J3 

romn(n,n2, 03)^2^3 



1+ 



X+ 



0(1) , fc>0 (B.26) 
1 log(-n)(F„n(n2, □3)^2^3) +0(-^. (B.27) 
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For all subcases in (B.5) one has 



TkmniO, ^2, □3)^2X3 



romn(n,n2, □3)^2X3 



1+ 



1+ 



= - '",I'ri'!:v'' f^4n.°3)A-.X3 + 0(-L) . (B.28) 



(k + m + n)\ 
(m + n)\ 



k>0 



log(-n)(F^„(n2, □3)^2^3) 



+ i^mn (02, □3)^2X3 



+ . (B.29) 



For the superpositions of F with 1/ □ one has 



1+ 



_ (-ir+"(fe-i)! 1 . . 



{k + m + ny. □ 
+ ^^^S#«'»»(°-°3)« + 0(-L) ,(B.30) 



A; > 1 



1 

□ 



rimn(n,n2,n3)J2J3 



2 (A; + m + n)! 

-T^^W H(^-n(n2,n3)J2J3) 

^^ni M-n)(GU°2,n3)J2J3) 

2 (m + n + 1)! L ^ ^ 

+ M^n(n2, ^3)^2^3 - Gmn{02, 03)^2^3 



(B.31) 



The senior terms of the latter expressions, proportional to !/□, cancel in the energy- 
momentum tensor [1]. 

Another useful relation: 



□ 



Ffc„^„(□,□2, 03)^2^3 



(B.32) 



□2^3 

is a consequence of Eq. (B.14). 

Finally, the relations (B.26)-(B.32) remain unchanged when multiplied by □2/^3 or 
n^/n2. Indeed, replacing the function T{q,x) with (□2/n3)jF(g, x) or (□3/n2)jF(g, x) 
doesn't change the behaviours (B.6) and (B.7). Using Eq. (A. 3), the multipher □2/^3 or 
□3/02 can be absorbed in any second-order or third-order form factor. 
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Appendix C. The operators log(— □) and !/□ for test 

functions singular at 



The behaviour (4.3) of the resolvent is vahd only for test functions having compact 
spatial supports [10]. Therefore, the behaviours at of all form factors, used in [2], 
are also valid only under conditions (1.43). For log(— □) this is the behaviour (3.25). In 

the general case, the support of the test function may conventionally be divided into a 
compact domain and asymptotic domain. The 1/r^ behaviour in (3.25) is a contribution 

of the compact domain. Any behaviour of log(— n)^ more singular than 1/r^ (call it 

1+ 

just singular) can only be a contribution of the asymptotic domain, i.e. of X 

Similarly, the regular behaviour of(l/n)X isl/r and is a contribution of the compact 
domain. A key to obtaining the contributions of the asymptotic domain is the fact that 
the null hyperplane reaches X+ at only one point of the celestial sphere [2]. Therefore, 
the singular contributions are always local in the angles although possibly nonlocal in 
time. To see why they may be nonlocal in time recall that, when a point tends to X"*", 
one generator of its past light cone merges with X+ entirely [2] . The retarded time ranges 
along this generator to — oo whereas the whole generator is labelled by a single value of 
the angles. 

Let C{x, x) be the retarded kernel of log(— □) , 

\og{-D)X{x) ^ J dxf/^jC{x,x)X . (C.l) 

By the argument above. 



log(-n)X(a:) 



Jdxg'/'jCix,x)lx\^^^^+0{^) , (C.2) 



X = {u,(f),r ^ oo) 

i.e. for obtaining the singular terms at X+, the test function can be taken at the angles 
of the observation point, = 0. The angle integrations in dx can then be done explicitly, 
and, as a result, the kernel becomes spherically symmetric. Thus, for obtaining the 
singular contributions at X+, it suffices to consider the spherically symmetric kernel. 
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The spherically symmetric kernel suffices for obtaining also the regular contributions 
provided that X is a scalar, and one needs only the integral of (C.2) over the 2-sphere. 
Indeed, to lowest order in the curvature, the scalar kernel of log(— □) can depend on the 
angles only through the arc length between the points (f) and 4> on the 2-sphere. Therefore, 

I (fS{(l)) (log(-n)x) = I dxg^/''C{x,x)(^j d^S{4>)X^ (C.3) 

where the angle integrations in dx concern only the kernel C{x,x) and convert it into a 
spherically symmetric kernel. 



One case of the singular behaviour considered below is where X 



0(l/r), and 



\og{-U)X 



1+ 



is needed up to the regular terms 0(l/r^). In this case Eq. (C.2) works. 



Another case is where X 



x+ 



is needed including the regular 



0(l/r2), and \og{-n)X 
terms This case is more difficult but is encountered only in ^^^(3) (Sec. 7) where 

the limitations implied in Eq. (C.3) are fuffiUed. Therefore, in both cases one may use 
the spherically symmetric kernel. 

Below, y is a point of the 2-dimensional Lorentzian section of a spherically symmetric 
spacetime, and Y{y) is a test function restricted to this section. The spherically symmetric 
retarded kernel of the operator log(— □) is of the form [8,5] 



r 

df Y 

path 1 



u _ 

log( -)Y{y) = - / df^Y 

oo 
r 



+ 



1 /" d 

- df log(m(r — f)) — (fY 
r J dr^ 



V 

1 /■ d - 

- I df log(m(f — r)) — [fY 
r J dr^ 



+ 

+ 2cy(y) 



path 2'' 



path 3'' 



(C.4) 



where r is the luminosity coordinate of the observation point y, and the integrations are 
along the null pathes 1,2,3 shown in Fig. 2. In (C.4), each of the pathes is parametrized 
with the luminosity coordinate f. The retarded time labelling the radial future light 
cones and normalized in (1.11) will be denoted u as above. In the coordinates y — {u, r), 
y — {u, f), and with the curvature neglected, path 1 is u + 2f — u, path 2 is u — u, and 
path 3 is u + 2f — u + 2r. 
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The spherically symmetric retarded kernel of the operator !/□ is of the form 

-^Yiy) = ^Jd'yg'/\y)fY (C.5) 
n 

where is the domain bounded by the pathes 1,2,3, and d?y g^^^{y) is the induced volume 
element. Hence, in the coordinates y — {u,r), 





d_ 1 
du 



1 If 1 f 1 f 

—Y(u,r) — — drrY +— / drrV — — / 

□ ^ ^ 2rJ pathl 2rJ path 2 2r J 



dr rY 



path 3 



(C.6) 



Denoting the contributions of the pathes 1,2,3 in (C.4) Vi, V2, V3, one has 
- log(--^)r(y) = VM + V2{y) + VM + 2cY{y) . 
The contribution of path 2 can be rewritten identically as follows: 



(C.7) 



V2{y) 



7 

/ 



drr 



r — r 



-Y 



M\ogmT)Y{y) + \ jdi log(l - 0/(rO (C.8) 
path 2 r" J 



with n arbitrary, and 



/(.) . 



,) ■ 



path 2^ 

If one chooses n equal to the power of decrease of Y at X+ 



x+ r" 



(C.9) 



(C.IO) 



the last integral in (C.8) will have a finite limit: 



di iog(i-0/K) 



-A{u) 



(C.ll) 



In this way one obtains for n = 1 

log(mr) 



n = 1 , V2{,y) 



x+ 



(C.12) 



(with no pure 1/r term), and for n = 2 



ri = 2 , ^2(2/) 



^log(mr) ,, , 1 /■ , , / s. d / ^- 
2 ^\ M m - — df log mf —(f^F 

^2 J gjj, V 



path 2 



logr 



(C.13) 
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As y — > X+, path 3 shifts entirely to X~^. Introducing the retarded time m as a 
parameter along path 3, one can easily calculate the limit 

u — u\ d 



- ( 

du log(m- 



— oo 



du 



) 



(C.14) 



Hence for the behaviour (C.IO) one obtains 



x+ 



— I dulog(m^-^—^)-^A(u) + o(^— 



(C.15) 



Finally, as y — > X"*", path 1 remains fixed. Therefore, its contribution is always regular: 

_ 



1+ 



— drrY 



path 1 



(C.16) 



The contributions of the pathes 1,2,3 in (C.6) are considered similarly. 
In the case n = 2 above, the total result is 



-log{-n)Y{y) 



2^(logr + c) + -^ + 0(^) 



B(u) 



' log r ^ 



where 



B{u) = Bi{u) + B2{u) , 

d 



BAu) = [drrY — [ dr logr —(r^Y ) 

I pathl I df^ path 2^ 



path 2^ 



iu) = Jdulog{'^)^Aiu) , 



(C.17) 

(C.18) 
(C.19) 
(C.20) 



and the next task is obtaining the behaviour of the coefficient B{u) as u ^ oo. 

The analysis of the behaviour of Bi (u) at late time essentially repeats the one in [2] . 
The dominant contribution to this behaviour comes from Y{y) at the limit y ^ which 
in the present case is the limit r ^ oo along the radial timelike geodesic that reaches the 
future infinity with the energy £^ = (1 — 7^)"^^^ per unit rest mass: 

1-7 



u 



7 



r , r — >• oo 



(C.21) 



The variables 7 and r may be used as coordinates of the point y : 

y(y)=y(7,r) . 
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(C.22) 



Then the definition of the hmit i"*" is 



Y 



F(7, r — > go) 



Of interest is the following behaviour of Y at i"*" (see Sec. 7): 

y = ^^^^—^Qii) , Q(i) 7^ 



(C.23) 



(C.24) 



where Q{'-f) is some regular function of 7. 

The limits i"^ and X+ are related [2]. For an analytic function, the sequence of limits 
i"*" and 7 — 1 coincides with the future of X"*". Hence, using (C.21), one obtains 



(F ) =(F 



(C.25) 



Therefore, the behaviour (C.24) implies a linear growth of the coefficient in (CIO) at late 
time: 

A{u) =2uQ{l) . (C.26) 

The late-time behaviour of Bi{u) in (C.19) is obtained by introducing 7 as an integra- 
tion variable in both integrals and restricting0 both integrations to the interval < 7 < 1. 
One obtains 

1 

(i7 



BAu] 



—u 



— u 



{1+1} 

I 



(1-7) 



-/ii(7,r 



7U 
1 + 7 

7M 



00 



00) 



where 



hi{'~f,r) = rY{'j,r) 



/ 9 7(1 -7) a N 2^ 



With the behaviour (C.24) of Y at i^ this yields the result 



(C.27) 

(C.28) 
(C.29) 



Bi(u] 



1 - (9 
u I d77Y^Q(7)-^ / rf7(log^) — (7'(l + 7)Q(7))+0(logn) . 
7 7 

(C.30) 



^"^The integration limits < 7 < 1 emerge after one restricts the support of Y to the interior of some 
future hght cone u = const, and the exterior of some tube f — const. The complementary portions of the 
support of Y contribute negligibly as u — > 00 [2] . 



57 



The integration by parts brings this expression to the final form 



^ 2 

-2M(log«-log2 + 2)g(l)+4M / rf7-^(g(7)-g(l))+0(logM) (C.31) 



in which the coefficient of the hnear growth is convergent. In this way the pole at 7 = 1 
is eliminated (cf. Eq. (1.21) and the discussion of this problem in [2]). 

The behaviour of B2{u) at late time is obtained by rewriting Eq. (C.20) identically as 
follows: 



— J 

B,{u)^ J dulog{^)-A{u) + {log^){A{u)-A{uo))+u J ci^ log(l - ^K) 

—00 uo/u 

(C.32) 

where uq < u, and 

g{u)^-^A{u) . (C.33) 

As u ^ 00, the first term in (C.32) is 0{logu), and the remaining terms are determined 
by the behaviour (C.26). In this way one obtains 



B2{u) 



2u{\og^-l)Qil) + 0{\ogu) 



(C.34) 



In the sum (C.18) the senior terms ulogu cancel, and the final result is 



B{u) 



u 



i 2 

-6Q(l)+4|d7^4-^(g(7)-Q(l))]+0(logt.) . (C.35) 



Taking into account Eqs. (C.2) and (C.3), one can summarize the calculations above 
as follows. For any function X{x) in four dimensions that behaves at like 

A{u,(f)) 



X 



n<2 



x+ 



one has 



X 

□ 



and 



-log( -)X 



x+ 2(2 - n) 

A(m,0) 



(C.36) 



(C.37) 



1+ 



log mr + 2c — log 2 — 1 + 

d 



dj 1-e 



+ 1/ cirlog(m(«-r))^A(r,0) + ^ (C.38) 
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where O 



— 0. For a function X(x) that behaves at X+ hke 



X 



A{u,<P) 



x+ 



one has 



and 



□ 



1 logr 



j drA(T,0) + O(-) 



(C.39) 



(C.40) 



\og{-U)X 



= 2 



(.0g. + c) + ^ + 0(!^^) (C.41) 



x+ 

with some coefficient B{u,(j)). If in the latter case the function X{x) is a scalar that 
behaves at i"*" like 

7(1-7^^ 



X 



(C.42) 



then 

d'^S{(j))B{u,(i)) 



u 



£S{<I>) [-6g(l,0)+4 J d^j^{Q{-f,<l>)-Qil, 



+ 0{\ogu) 



(C.43) 



Note that the term with logr in (C.41) is doubled as compared to (C.38). 
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Figure captions 



Fig.l. The diagram for the contribution (B.l) to the expectation- value current at point x. 
The argument □ of the form factor T corresponds to the external line. 

Fig. 2. Penrose diagram for the Lorentzian section of a spherically symmetric spacetime. 
The timelike line r = is the central geodesic. The union of pathes 1 and 2, and 
path 3 are the two radial light rays that come to the 2-dimensional observation 
point y. fl is the domain bounded by the pathes 1,2,3. 
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Fig.l 
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